Localization of bosonic atoms by fermionic impurities in a 3d optical lattice 
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We observe a localized phase of ultracold bosonic quantum gases in a 3-diniensional optical lattice 
induced by a small contribution of fermionic atoms acting as impurities in a Fermi-Bose quantum 
gas mixture. In particular we study the dependence of this transition on the fermionic '"^K impurity 
concentration by a comparison to the corresponding superfluid to Mott insulator transition in a pure 
bosonic *^Rb gas and find a significant shift in the transition parameter. The observed shift is larger 
than expected based on a mean-field argument, which is a strong indication that disorder-related 
effects play a significant role. 
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Ultracold atomic gases have already shown to be ver- 
satile model systems for quantum phenomena in many 
areas of physics. Prominent examples are the recent 
demonstration of the superfluid to Mott insulator tran- 
sition [l| and fermionic superfluidity in the BEC-BCS 
crossover regime Q, |^ |j, |M ly| ■ These phenomena of so- 
called "strongly correlated physics" originate from the 
rich field of transport and localization in condensed mat- 
ter systems with high relevance for technical applications. 
Realizations of such systems with ultracold atoms in op- 
tical lattices are a breakthrough as they offer a "defect- 
free" system as well as in situ control over tunneling 
and interaction parameters, allowing a direct quantita- 
tive comparison to theory. However, in contrast to single- 
component quantum gases, real condensed matter sys- 
tems are not defect-free and impurities play a crucial role, 
namely causing important localized phases connected to 
charge-density waves, disorder or percolation. Thus, the 
realization of impurity physics with quantum gases is of 
high relevance. Fermi-Bose mixtures in 3-dimensional op- 
tical lattices provide a versatile model system for many 
disorder- and impurity-related effects. 

Interest in this novel strongly correlated mixed statis- 
tics system has been driven by investigations of new 
Fermi-Bose correlations and their consequences |3, |sl Isl 
|liJ,|ll|,|l3,|lJ|. In the extreme case of pairing of fermions 
with one or more bosons a whole zoo of new quantum 
phases of these "composite fermions" has beenpredictcd 
jl4j , and connections to superconductivity 0, lq| but also 
"exotic" systems like neutron stars and quark-gluon plas- 
mas [l^ suggest fundamental relevance for a broad range 
of phenomena. Even before such atom pairs form, Fermi- 
Bose correlations are predicted to become manifest in 
the formation of charge density waves [l^ llTl Il8| or 
supersolids [li|, but also in polaron-related physics of 
fermions dressed with bosons [la| or percolation [ijj. 
Such disorder-driven localization phenomena are present 
and intensely discussed in many systems, from Anderson 
localization in condensed matter and light localization in 
random media to seismic wave dynamics and even dis- 
turbances in traffic flow. 

In this letter, we present first studies on Fermi-Bose 




FIG. 1: Schematic localization scenarios, a. Pure bosonic 
superfiuid in an optical lattice, b. Shift of the effective po- 
tential depth due to fermionic impurities, c. Localization 
by interfering paths of the bosonic wavefunction scattered by 
randomly distributed fermionic impurities, d. Localization 
due to percolation. A random fermion distribution hampers 
the establishment of a coherent connection and causes the 
localization of bosonic ensembles in superfluid "islands", e. 
Mott insulator transition induced by a uniform distribution 
of attractive fermionic impurities, resulting in an effectively 
deeper lattice potential for the bosons. 



mixtures confined in 3-dimcnsional optical lattices where 
fermionic atoms act as impurities. We observe the on- 
set of a localized phase in the bosonic ^^Rb component 
induced by the fermionic *'^K impurity atoms. We find 
a fermion concentration dependent shift with respect to 
the quantum critical point of a pure bosonic ensemble 
which is incompatible with simple mean-field arguments. 
This indicates that the observed transition is driven by 
disorder-related many-body effects. The measurements 
presented in this letter are the first to combine three- 



dimensional crystals made of light with quantum parti- 
cles of mixed statistics. On the one hand, this novel sys- 
tem is a realization of the Bose- Fermi Hubbard Hamilto- 
nian [lll|20j; on the other hand, it offers further potential 
for controlled investigations of impurity-induced phases 
as well as studies on pairing between bosons and fermions 
leading to composite fermions ll4L[l5l |21| and BCS-like 
boson-assisted fermion pairing p, la] 

We first consider intuitive localization scenarios for 
Fermi-Bose mixtures confined in 3-dimensional optical 
lattices, starting with a pure bosonic superfluid (Fig.ni). 
Adding fermionic impurities and considering the attrac- 
tive *°K — *^Rb Fermi-Bose interaction energy as an 
additional potential for the bosons, the "defects" caused 
by the fermionic impurities can be described by a local 
change of the effective optical lattice depth for the bosons 
due to the interparticle interaction (Fig. ^). If the en- 
ergy level shift caused by the interaction energy is large 
enough, the superfluid bosonic wavefunction will not ex- 
tend into this defect region, but will be scattered by the 
impurity. If scattering becomes frequent, interference ef- 
fects along a closed scattering path arc predicted to sup- 
press transport and lead to a localization scenario similar 
to Anderson localization |23| (Fig.^^). A further increase 
in the impurity density will lead to the formation of "for- 
bidden walls" . Once the walls in this quantum percola- 
tion scenario lead to a sufficiently complicated labyrinth 
like structure for the bosonic wavefunction, a single co- 
herent superfluid phase can no longer be sustained and 
several separated domains will be formed (Fig.njl). For 
the maximum filling of one fermion per lattice site the 
localized phase becomes comparable to a pure bosonic 
Mott-insulator but now the transition is shifted by the 
interaction energy with the fermionic impurity (Fig.^). 
Access to these various regimes in the experiment can be 
controlled by varying both crystal depth and fermionic 
impurity concentration. 

Our experimental procedure for creating degenerate 
mixtures has been described previously |23l |2J| . In brief, 
we create Fermi-Bose mixtures by first collecting and pre- 
cooling bosonic ®^Rb and fermionic ''"K atoms using a 
combined 2d-3d MOT laser cooling system. After trans- 
fer to a loffe-Pritchard type magnetic trapping potential 
we perform forced RF evaporative cooling of the rubid- 
ium component in the |2,2) state (which in turn sym- 
pathetically cools the potassium atoms in the |9/2, 9/2) 
state) until the mixture is close to quantum degeneracy. 
At this point we add a dipole trapping potential in or- 
der to increase the axial trapping frequency from 11 to 
wSO Hz (this and all subsequent trap frequencies cor- 
respond to the ^'^Rb component) and at the same time 
reduce the radial trapping frequencies from 250 Hz to 
150 Hz, resulting in a significant decrease of the trap as- 
pect ratio. Having created a quantum degenerate Fermi- 
Bose mixture by a last evaporation step in this combined 
trap, either this mixture is used for further experiments 




FIG. 2: Time of flight absorption images of the bosonic com- 
ponent 15 ms after switching oH' the lattice and trap poten- 
tials. (The colors and pseudo-3d representation encode the 
particle column density integrated along the imaging direc- 
tion), a. Pure bosonic ensemble for three different lattice 
depths, b. Fermi-Bose mixture for the same lattice depths 
and 20% fermionic impurity concentration. The reduction of 
interference contrast visible in the images shows the onset 
of localization, which is shifted to lower lattice depth in the 
presence of fermionic impurities. 



or the fermions are removed with a short resonant light 
pulse, leaving a pure bosonic sample as reference. The 
ultracold mixture of about 10^ bosonic *^Rb atoms with 
no discernible thermal cloud contains a variable 5-20% 
fermionic '"'K atom impurity component. 

The optical lattice potential is created by three orthog- 
onal and retroreflected laser beams (waists 82/92 /im in 
the radial directions and 55 fim in the axial direction 
with respect to the original trap) with mutually orthog- 
onal polarizations and a relative frequency difference of 
about lOMHz. The lattice is operated at a wavelength 
of 1030 nm, creatin g a simple cubic crystal with a lattice 
constant of 515 nm J25|. Due to the additional harmonic 
confinement by the magnetic trap and the Gaussian lat- 
tice laser profiles, the mixture occupies a few ten thou- 
sand lattice sites with an occupation rising from in the 
outer regions to 1 fermion and >5 bosons per site at the 
center. Note that we have chosen parameters such that 
the fermionic impurities always stay within the bosonic 
cloud. After sudden switch-off of the lattice potential 
and a period of free expansion of 15-20 ms we record the 
interference pattern which builds up in the density dis- 
tribution of the bosonic component [26j . In order to map 
the localization transition we repeat this procedure for 
several different final lattice depths. 

Fig. [21 shows the evolution of the boson interference 
pattern for about 20% impurity concentration (bottom 
row) in comparison to a pure bosonic sample created 
under the same experimental conditions but with a re- 
moval of the fermionic atoms just before the lattice ramp- 
up sequence (top row). The lattice depths are given 
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FIG. 3: Visibility V (a) and central peak width (b) of the 
bosonic interference pattern for different fermionic impurity 
concentrations. 



in units of the recoil energy for the ®^Rb component 
E^^ = {h^k^)/{8TT^mmy) « h ■ 2.14kHz where k is the 
lattice wavenumber. The loss of interference contrast 
marking the breakdown of long range order is clearly vis- 
ible in both cases. In case of the pure bosonic gas, the 
loss of coherence is due to the well-known superfluid to 
Mott-insulator phase transition [l|, |23] which occurs as a 
result of competition between the minimization of kinetic 
energy, parametrized by the tunneling matrix element J 
which tends to delocalize the atomic wavefunction over 
the crystal and the minimization of interaction energy U 
(Fig. |2Ji). As can be clearly seen from Fig. 13), the pres- 
ence of fermionic impurities induces a loss of coherence 
at much lower lattice depths than for a pure BEC. 

In the following, we will study this impurity-induced 
phenomenon quantitatively. Taking the interference con- 
trast in the above pictures as a measure of coherence and 
putting it into mathematical terms, one can define the 
visibility of the matter wave interference pattern as [23 

V = iNpeaks-Nint)/{Npeaks+Nint)- Heie Npeaks deuOtCS 

the sum of the number of atoms within the first order in- 
terference maxima, while Ni„t is the sum of the number 



of atoms in equivalent areas at intermediate positions 
between the maxima. In addition, we extract comple- 
mentary information from the width of the central inter- 
ference peak as a measure of the correlation length of the 
bosonic system [l|,|23. Fig. |31 shows sample data compar- 
ing the behavior of ensembles with 3% and 7% impurity 
concentrations to a pure bosonic ensemble. The onset of 
localization marked by a loss of interference contrast and 
a decrease of correlation length is clearly visible in the 
visibility curve as well as in the width data. Fig.Ocleaiiy 
shows that the onset of the localized phase is shifted by 
an amount depending on the impurity concentration. Al- 
ready for 3% impurity concentration, the visibility data 
shows a significant shift. The corresponding width data, 
however, essentially exhibits a steeper slope in the local- 
ized phase with no pronounced shift as compared to the 
pure bosonic system. This indicates a qualitative differ- 
ence in the type of transition for the pure bosonic system 
and the system with fermionic impurities. 

In order to assess the size of the shift, one can calculate 
the value of the on-site Fermi-Bose interaction matrix 
element 
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The value on the right-hand side is calculated in the 
tight binding limit at the localization transition point 
for our lattice parameters, based on a scattering length 
of ofb = — 205(7)ao J3fl |. This matrix element charac- 
terizes the additional potential felt by bosonic atoms due 
to the presence of a fermionic impurity atom in a single 
potential well. It thus corresponds to the expected shift 
of the localization transition for unity fermionic filling in 
the simple mean field picture shown in Fig. ^. 

We use two different methods to extract the shift of the 
localization transition from our experimental data. For 
the visibility data, we use a phenomenological fit function 



V(,s) 



V(0) 



1 + Exp(q; • (s - Scrit)) 



where V(0) is the visibility for shallow lattices, s is the 
lattice depth in units of E^^ , Scrit is a measure for the 
onset of the localization transition and a is an additional 
fit parameter. For data as in Fig. |31 we compare Scrit 
for a pure bosonic sample recorded under the same ex- 
perimental conditions to the impurity induced transition 
point and extract the shift of the transition. Resulting 
shifts are plotted as a function of impurity concentration 
in Fig. 21 The corresponding width data is analyzed by 
extracting the intersection point of two linear fits to the 
descending and ascending branches of the data. As can 
be seen from Fig. ^ an increasing impurity concentra- 
tion leads to a considerable shift of the localization tran- 
sition, significantly larger than the on-site Fermi-Bose 
interaction energy of C/fb- For 20% impurity concentra- 
tion, the shift is on the order of hE^"" , about 5 times 
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FIG. 4: Observed shift of the localization transition as a func- 
tion of fermionic impurity concentration 



t/pB • Our data is therefore clearly incompatible with the 
above mean-field argument even for unity fermionic fill- 
ing |31|, and indicates additional effects like disorder. At 
the same time, for smaller impurity concentrations, the 
shift extracted from the width is always smaller than the 
visibility data. Thus our data suggest that global phase 
coherence is reduced with fewer impurities than the cor- 
relation length, also indicating the presence of disorder 
effects. 

We have carefully checked temperature and ramp 
speed dependence as well as heating and loss effects 
and found that despite some losses in fermion number 
compatible with three body recombination |23l | and mild 
heating, the impurity effects caused by the fermions are 
clearly identifiable. Furthermore the localization transi- 
tion is reversible, i.e. the visibility is nearly fully recov- 
ered, when the lattice depth is ramped back down, which 
shows that non-adiabatic effects play a minor role. 

Fermi-Bose mixtures in a three-dimensional optical 
crystal as presented in this letter open many new op- 
portunities for the experimental investigation and under- 
standing of localization phenomena on various grounds. 
Future studies, e.g. using second order correlations or 
Bragg spectroscopy, will concentrate on the further iso- 
lation and identification of the underlying phases and as- 
sociated phase transitions in this rich system. 
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